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We consider massless higher spin gauge theories with both electric and magnetic sources, with 
a special emphasis on the spin two case. We write the equations of motion at the linear level 
(with conserved external sources) and introduce Dirac strings so as to derive the equations from a 
variational principle. We then derive a quantization condition that generalizes the familiar Dirac 
quantization condition, and which involves the conserved charges associated with the asymptotic 
symmetries for higher spins. Next we discuss briefly how the result extends to the non linear theory. 
This is done in the context of gravitation, where the Taub-NUT solution provides the exact solution 
of the field equations with both types of sources. We rederive, in analogy with electromagnetism, 
the quantization condition from the quantization of the angular momentum. We also observe that 
the Taub-NUT metric is asymptotically flat at spatial infinity in the sense of Regge and Teitelboim 
(including their parity conditions). It follows, in particular, that one can consistently consider in 
the variational principle configurations with different electric and magnetic masses. 

I. INTRODUCTION 

The idea of electric-magnetic duality has resiliently resisted developments and appears to be here to stay. Originally 
explored by Dirac for electromagnetism [1, 2], it was later analyzed in the context of non-abelian gauge theories in 
[3-5]. It was more recently generalized to extended objects and p-form gauge fields in [6, 7]. A fascinating implication 
of electric-magnetic duality is the charge quantization condition. This condition is antisymmetric for p-dyons of even 
spatial dimension p, and symmetric for odd p [8-10]. 

Fields with spin higher than that of the graviton were discussed some time ago by Fronsdal [11] who gave the action 
integral for the free theories (see also [12]). Interest in these higher spin fields has been renewed since they appear in 
the spectrum of string theory. In particular, duality for free higher spins has recently been the focus of many papers 
[13-18]. Conserved, external, electric-type sources can easily be coupled to the higher spin fields but the problem 
of constructing consistent self-interactions is still incompletely understood beyond spin 2 (for instance, the action 
principle is not known even though non linear equations have been constructed in the remarkable work [19]). It has 
been argued that self-interactions can be brought in but only among an infinite tower of fields [19]. 

Since duality can be defined for higher spins and since electric sources can be included, one might wonder whether 
magnetic sources can be considered as well. Our paper solves positively this question for all spins at the linearized 
level and provides additional insight in the full non linear theory for spin 2. 

We show that conserved external sources of both types can be coupled to any given higher (integer) spin field 
within the context of the linear theory. The presence of magnetic sources requires the introduction of Dirac strings, 
as in the spin 1 case. To preserve manifest covariance, the location of the string must be left arbitrary and is, in 
fact, unobservable classically. The requirement that the Dirac string is unobservable quantum- mechanically forces a 
quantization condition of the form 

lQ 1r . Vl WF-^We2. (LI) 

Here, the symmetric tensor P 71 ' 7a ~ 1 (m) is the conserved electric charge associated with the asymptotic symmetries 
of the spin s field, while Q 7l ... 7a „ 1 (v) is the corresponding "topological" magnetic charge. For s = 1, the asymptotic 
symmetries are internal symmetries and, actually, just constant phase transformations. The conserved charge P is the 
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electric charge q while Q is the magnetic charge g, yielding the familiar Dirac quantization condition for the product of 
electric and magnetic charges. For s = 2 the conserved charges have a spacetime index and the quantization condition 
reads (after rescaling the conserved quantities so they have dimensions of mass) 

EZ. (1.2) 

h 

The quantity P 1 is the "electric" 4-momentum associated with constant linearized diffcomorphisms (translations) 
while Q 7 is the corresponding magnetic 4-momentum. For a point particle source, P 7 = Mm 7 where M is the 
"electric" mass and u 7 the 4-velocity of the electric source. Similarly, Q 1 = Nv 7 where N is the "magnetic" mass 
and i> 7 the 4-velocity of the magnetic source. 

All this is just a generalization of the familiar spin 1 case, although the explicit introduction of the Dirac string 
is more intricate for higher spins because the gauge invariance is then more delicate to control. Indeed, there is no 
gauge invariant object that involves first derivatives of the fields only (s > 1). Hence, the Lagrangian is not strictly 
gauge invariant, contrary to what happens for electromagnetism, but is gauge invariant only up to a total derivative. 

A serious limitation of the linear theory for s > 1 is that the sources must follow straight lines. This follows from the 
strict conservation laws implied by the field equations, which are much more stringent for s > 1 than they are for s = 
1. Thus the sources must be treated as externally given and cannot be freely varied in the variational principle. One 
cannot study the back reaction of the spin s field on the sources without introducing self-interactions. This problem 
occurs already for the spin 2 case and has nothing to do with the introduction of magnetic sources. 

We have not investigated the backreaction problem for general spins s > 2 since the non linear theory is still a 
subject of investigation even in the absence of sources. We have focused instead on the spin 2 case, for which the non 
linear theory is given by the Einstein theory of gravity. The remarkable Taub-NUT solution [20] , which represents the 
vacuum exterior field of a gravitational dyon, indicates that Einstein's theory can support both electric and magnetic 
masses. The remaining part of our paper is devoted to a further discussion of the Taub-NUT solution, including an 
analysis of its asymptotic properties. We provide new insight on the quantization condition (1.2) through the analysis 
of the angular momentum stored in the gravitational field of a pair of electric and magnetic masses placed at different 
locations. 

Our paper is organized as follows. In section II, we consider in detail the spin 2 case with point particle electric and 
magnetic sources. We introduce Dirac strings and derive the quantization condition. We then extend the formalism 
to higher spins (section III), again with point particle sources. In section IV we consider the Taub-NUT solution and 
discuss the quantization condition in the non linear context. After a brief concluding section, the article ends with 
three appendices where technical aspects for spins greater than 2 that were not included in section III are dealt with 
(Appendices A, B and C). Finally some aspects of the angular momentum of an electric test particle in the field of a 
monopole, which turn out to be relevant to the discussion of the Taub-NUT metric are covered in Appendix D. 



II. LINEARIZED GRAVITY WITH ELECTRIC AND MAGNETIC MASSES 



A. Electric and magnetic sources 

The equations of motion for linearized gravity coupled to both electric and magnetic sources are naturally written 
in terms of the linearized Ricmann tensor R a fjXfn hereafter just called "Ricmann tensor" for simplicity. This is the 
physical, gauge-invariant, object, analogous to the field strength of electromagnetism. How to introduce the 
"potential", i.e., the symmetric spin 2 field = h ufl will be discussed below. The dual to the Ricmann tensor is 
defined as 

We denote the "electric" energy- momentum tensor by T^ v and the "magnetic" energy- momentum tensor by 9 Miy . 
These are both symmetric and conserved, T^ v = T v », 9^ = 9^, T» v v = 0, ®» v v = 0. It is also useful to define 
f ^ _ rp^iu _ i Tj q^u = e ^ _ i q where T and are the traces . We assume that T* v and 9^ have the 
units of an energy density. We set c — 1 but keep G. 

The form of the equations in the presence of both types of sources is fixed by: (i) requiring duality invariance with 
respect to the 50(2)-rotations of the curvatures and the sources [17], 

R'apXfj, = cos a R a p\ l j, + sin a Sapxfi, S' a/3X ^ = - sin a i? Q/3 A M + cos a S a pxn, 
T' a[j = cos a T a fj + sin a 9 a/ g, 9' a/3 = - sin a T a p + cos a 9 Q( a, 
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and, (ii) using the known form of the equations in the presence of electric masses only. One finds explicitly the 
following: 

• The Riemann tensor is antisymmetric in the first two indices and the last two indices, but in general is not 
symmetric for the exchange of the pairs, i.e., R a px^ = -RpaX^, R a fix^ = -Rap^x with R afi x p ^ R\^ a [) (in the 
presence of magnetic sources). 

• In the presence of magnetic sources the cyclic identity is [33] 

Ra/3\fJ. + R/3\an + R\af3n = 87rG e Q/3A „ 0"^ . (H.l) 
This enables one to relate R a px^ to Rx^ap through 

Ra/3-fS - RySafJ = 47rG (e Q /3 7 A@ A ,5 ~ e Q /3<5A@ A 7 + £/3 7 <5aB A q - e al S\0 X p) ■ (II. 2) 

It follows that the Ricci tensor is symmetric, R\^ = R^\. The Einstein tensor = R\^ — (l/2)rj\ ll R is then 
also symmetric. 

• The Bianchi identity is 

d e R a f3j8 + d a Rp e ^s + dpR^s = 87rG e eQ(3 p(9 7 9 p (5 - d$<d p ^) . (II. 3) 

Although there is now a right-hand side in the Bianchi identity, the contracted Bianchi identities are easily 
verified to be unaffected and still read 

G V AJ =0. (II.4) 

• The Einstein equations are 

G x » = 8irGT Xtl , (II.5) 

or equivalently, R Xfl = 8ttGT Xi1 , and force exact conservation of the sources because of the contracted Bianchi 
identity, as in the absence of magnetic mass. 

The equations are completely symmetric under duality. Indeed, one easily checks that one gets the same equations 
for the dual curvature S a px^ with the roles of the electric and magnetic energy-momentum tensors exchanged. In the 
course of the verification of this property, the equation 

d^R^s = SttG (d^f p s - d s f pl ) , 

which follows from (II. 2), (II. 3) and the conservation of are useful. Furthermore, in the absence of magnetic 
sources, one recovers the equations of the standard linearized Einstein theory since the cyclic and Bianchi identities 
have no source term in their right hand sides. 

The formalism can be extended to include a cosmological constant A. The relevant curvature is then the MacDowell- 
Mansouri curvature [21] linearized around (anti) de Sitter space [18]. In terms of this tensor, the equations (H.l), 
(II. 3) and (II. 5) take the same form, with ordinary derivatives replaced by covariant derivatives with respect to the 
(anti) de Sitter background. 

B. Decomposition of the Riemann tensor - Spin 2 field 

We exhibit a variational principle from which the equations of motion follow. To that end, we first need to indicate 
how to introduce the spin 2 field h pv . 

Because there are right-hand sides in the cyclic and Bianchi identities, the Riemann tensor is not directly derived 
from a potential h^. To introduce h^, we split Rx^ap into a part that obeys the cyclic and Bianchi identities and a 
part that is fixed by the magnetic energy-momentum tensor. Let <& Q/3 7 be such that 

^$"^ = 167^6^, $^ 7 = -$^ 7 . (II.6) 

We shall construct <& Q/3 7 in terms of 0^ 7 and Dirac strings below. We set 

Rx^p = r x »a(} + \ ex wa (5 Q $ P % - d^ p \) , (II.7) 
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with 

$> p \ = $ pff Q + i (<5 P Q $ CT - <T Q $ P ) , $ p = $ pff CT . 

Using d a $ Q/3 7 = IGttG©^ — <9 7 i >/3 , = — \^ , one easily verifies that the cyclic and Bianchi identities take the 
standard form when written in terms of r Q/ 3A M , namely, 

rapxn + rp\an + rxafj^ = 0, d e r a i3 7 s + d a rp ej s + dpr^s = 0. 

Hence, there exists a symmetric tensor such that r a [3Xfi = — 2d[ph a ][\. ^ . 
If one sets y AAI 7 — e App<T <9 p /i CT7 = — y pA 7 , one may rewrite the curvature as 

Rx^p = \ ex W a (d a Y p % - d Y^ a ) , (II.8) 

with 

y p % = y p % + = - Ya % Y pa a = Y pa a + \{5 p cX a -KY p ), y p ^y pa a , (n.9) 

(note that yP° a = yP a a and that d p yP° a = 0). 



C. Dirac string 

We consider point particle sources. The particles must be forced to follow straight lines because of the conservation 
equations T plJ v = and Q^ v v = 0. If w p is the 4- velocity of the electric source and w M the 4- velocity of the magnetic 
source, one has 

T^ = Mu v j dX5^(x- z(A))i p , e"„ = Nv u j d\5 iA Hx - z{\))^ , (11.10) 

where z^(X) and z^(X)) are the worldlines of the electric and magnetic sources respectively, e.g. u p = dz^/ds. 
Performing the integral, one finds 



TV = ^fs^(x-z(x )), &»> = ^>>(x - i(x )) . 



The tensor $ a/3 7 introduced in (II. 6) can be constructed a la Dirac [2], by attaching a Dirac string y p (A, a) to the 
magnetic source, y^(X,0) = z p (A). One has explicitly 

= 16ttGNv 7 J dXdoiy^y 13 -y a y' p )5 iA Hx-y(\o)), (11.11) 

where 

V ~ d\ ' y dc 

One verifies exactly as for electromagnetism that the divergence of $ Q/3 7 is equal to the magnetic energy-momentum 
tensor (up to the factor \QttG). What plays the role of the magnetic charge g in electromagnetism is now the conserved 
product Nv p of the magnetic mass of the source by its 4-velocity. This is the magnetic 4-momentum. 



D. Variational principle 

1. Action 



When the curvature is expressed in terms of h pv as in (II. 8), the expressions (II. 1) and (II. 3) are identically fulfilled 
and the relations (II. 5) become equations of motion for h pv . These equations can be derived from a variational 
principle which we now describe. 
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The action that yields (II. 5) is 

S[h^(x),y»(\, a)} = ^ J 1 (F a/37 y Q ^ - Y a Y a ) d 4 x + ±J h^T^d 4 x . (11.12) 

One varies the fields h pv and the coordinates of the string (with the condition that it remains attached to the 
magnetic source), but not the trajectories of the sources, which are fixed because of the conservation laws dpT^" = 
and d^Q^ = 0. This is a well known limitation of the linearized theory, present already in the pure electric case. To 
treat the sources as dynamical, one needs to go beyond the linear theory. 
If there is no magnetic source, the first term in the action reduces to 

s PF = J \ (-d x h af3 d x h al3 + 2d x h Xa d^ a - 2d x hd^ x + d x hd x h) , 

which is the Pauli-Fierz action, whose variation with respect to h a p gives — 16 ] rG times the linearized Einstein tensor 
G a/3 . It is straightforward to verify that the variation of the first term in the action with respect to h a p still gives 
~I(37rG times the linearized Einstein tensor G a/3 with correct <&^ v a contributions even in the presence of magnetic 
sources. So, the equations of motion that follow from (11.12) when one varies the gravitational field are the Einstein 
equations (II. 5). 

Extremization with respect to the string coordinates docs not bring in new conditions provided that the Dirac 
string does not go through an electric source (Dirac veto). 

The action (11.12) was obtained by using the analysis of source-free linearized gravity in terms of two independent 
fields given in [16], which enables one to go from the electric to the magnetic formulations and vice- versa, by elimination 
of magnetic or electric variables. As one knows how to introduce electric sources in the electric formulation, through 
standard minimal coupling, one can find how these sources appear in the magnetic formulation by eliminating the 
electric variables and keeping the magnetic potentials. So, one can determine how to introduce electric poles in the 
magnetic formulation, or, what is equivalent, magnetic poles in the electric formulation. 

2. Gauge invariances 

Diffeomorphism invariance 

The action (11.12) is invariant under linearized diffeomorphisms and under displacements of the Dirac string (ac- 
companied by appropriate transformations of the spin 2 field). The easiest way to show this is to observe that the 
first term in the action (11.12) is invariant if one shifts Y tlv a according to 

Y^ a -» Y^ a + 5^d p z v P - 5»d p z»P + d a z^ , (11.13) 

where z^ v = —z v ^ is arbitrary. This is most directly verified by noting that under (11.13), the tensor Y tiv a defined in 
(II. 9) transforms simply as 

Y^ a ^Y^ a +d a z^ (11.14) 

and this leaves invariant the first term in (11.12) up to a total derivative. Note that the Riemann tensor (II. 8) is 
strictly invariant. The transformation (11.13) can be conveniently rewritten as 

Y^ a - Y^ a + e^d p a aa , (11.15) 

where a aa = -a aa is given by a aa = ^e a af3jZ^. 
A (linearized) diffeomorphism 

hp„ -» h pv + d^ v + d,,^ (11.16) 

(with the string coordinates unaffected) modifies Y^ v a as in (11.15) with a aa — <9 Q £ CT — <9 CT £ Q (note that the term 
d a S, a in a aa does not contribute because d[ p d a ]^ a = 0). Hence, the first term in the action (11.12) is invariant under 
diffeomorphisms. The minimal coupling term is also invariant because the energy- momentum tensor is conserved. It 
follows that the complete action (11.12) is invariant under diffeomorphisms. 
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Displacements of the Dirac string 

An arbitrary displacement of the Dirac string, 

y a (T,a)^y a (T,a)+dy a (r,a) (11.17) 

also modifies Y p,v a as in (11.15) provided one transforms simultaneously the spin 2 field h^ v appropriately. Indeed, 
under the displacement (11.17) of the Dirac string, the quantity &^" a changes as — > <fr pv a + k plJ a where k^ v a 
can be computed from Sy a (r 7 a) through (11.11) and has support on the old and new string locations. Its explicit 
expression will not be needed. What will be needed is that it fulfills 

d^\ = 0, (11.18) 

because the magnetic energy-momentum tensor is not modified under a displacement of the Dirac string. The field 
Y pl/ a changes then as 

Y^ a -» Y pv a + e^dpShva + hT a (11.19) 

where 8h aa is the sought after variation of h aa . By using (11.18), one may rewrite the last term in (11.19) as d p t pup a 

for some t tlup a — t^ p ^ a . Again, we shall not need an explicit expression for t plJp al but only the fact that because 
k^ v a has support on the string locations, which do not go through the electric sources (Dirac veto), one may choose 
t tiVp a to vanish on the electric sources as well. In fact, one may take t pup a to be non- vanishing only on a membrane 
supported by the string. Decomposing V lvp a as t pil/p a = e pvpa [s aa + a aa ), s aa = S( CTa ), a aa = a[ aa ] and taking h aa to 
transform as h aa —> h aa — s aa one sees from (11.19) that the variation of Y^ v a takes indeed the form (11.15). Thus the 
first term in the action is invariant. The minimal coupling term is also invariant because the support of the variation 
of the spin 2 field does not contain the electric worldlines. 

One can also observe that the variation 8r a p pa vanishes outside the original and displaced string locations. This 
implies 6h a p — d a ^p + d^ a except on the location of both strings, where £ Q induces a delta function contribution on 
the string ("singular gauge transformation"). The explicit expressions will not be given here. 

Identities 

The identities which follow from the invariance (11.13), or (11.15), of the first term 



in the action may be written as 



They imply that 



9 P I I ^ - 9 P (4^c ) ^ ■ (H.20) 



P \dY a(3 J P \dY^ 



' = = J^-W (n.21) 



16ttG Sh al3 \dY' 

^ (IL22) 

from which the contracted Bianchi identities are easily seen to indeed hold. 

The expression (II. 8) of the Ricmann tensor in terms of Y p,v a makes it clear that it is invariant under (II. 14) and 
thus, invariant under both diffeomorphisms and displacements of the Dirac string. 



E. Quantization condition 



Because of the gauge invariances just described, the Dirac string is classically unobservable. In the Hamiltonian 
formalism, this translates itself into the existence of first-class constraints expressing the momenta conjugate to the 
string coordinates in terms of the remaining variables. Demanding that the string remains unobservable in the 
quantum theory imposes a quantization condition on the charges, which we now derive. The argument follows closely 
that of Dirac in the electromagnetic case [2] . 
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Working for simplicity in the gauge y° = X (which eliminates y° as an independent variable), one finds the constraints 

dC 

^^-^GNy'-v,— . (11.23) 

The right hand side of (11.23) generates the change of the gravitational field that accompanies the displacement of 
the Dirac string. 

In the quantum theory, the wave functional ip must therefore fulfill 

ft Sip „„ „ » T /„ dC 

-iSy^) = - 327rGNy ^F^" 

We integrate this equation as in [2], along a path in the configuration space of the string that encloses an electric 
source. One finds that the variation of the phase of the wave functional is given by 



WirGNv-y 

A* = - 

n 



J ^(y m y' n -y n y ,m )d*d\, (n.24) 



where the integral is taken on the two-dimensional surface enclosing the electric source. Using the Gauss theorem, 
this can be converted to a volume integral, 



A* = - 1 / d x e mnp d p t^tt — 



Because e mnp d p ^ a y£„ = 3^-, the variation of the phase becomes, upon use of the constraint (initial value) Einstein 
equations G° 7 = 87rGT° 7 , 

A* = — - / d A xT Ul = — . 

h J h 

For the wave functional to be single-valued, this should be a multiple of 27r. This yields the quantization condition 

ANMGv^vTi 

— y - — =n, neZ. (11.25) 

Introducing the conserved charges P 7 , Q 1 associated with the spin 2 theory (electric and magnetic 4- momentum), 
this can be rewritten as 



4GP 7 Q 
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7 ^ eZ. (11.26) 



h 

It is to be stressed that the quantization condition is not a condition on the electric and magnetic masses, but rather, 
on the electric and magnetic 4-momenta. In the rest frame of the magnetic source, the quantization condition becomes 

4GEN 

^j— GZ, (11.27) 

where E is the (electric) energy of the electric mass. Thus, it is the energy which is quantized, not the mass. 

We have taken above a pure electric source and a pure magnetic pole. We could have taken dyons, one with charges 
(P 7 ,Q 7 ), the other with charges (P 7 ,Q 7 ). Then the quantization condition reads 

AG (P 7 7 - P 7 Q 7 ) _ 4G , e ab Q°Q fc7 

h = h GZ ' (IL28) 

since the sources are pointlike (0-dyons). Here = (P 7 ,Q 7 ), a, b = 1,2 and e a fc is the S'0(2)-invariant Levi-Civita 
tensor in the 2-dimcnsional space of the charges. 
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F. One-particle solutions 

1. Electric mass 

We consider a point particle electric mass at rest at the origin of the coordinate system. The only non-vanishing 
component of its electric energy momentum tensor is T 00 (x°,x) = MS^(x) while 8^" vanishes. There is no Dirac 
string since there is no magnetic mass. The metric generated by this source is static. The linearized Einstein equations 
are well known to imply in that case the linearized Schwarzschild solution, namely in polar coordinates 

/loo = = n rri other components vanish, 



or in Cartesian coordinates 



2GM , 2GM 
hoc = , hij — — g — XiXj, other components vanish. 



Indeed, one then finds 



/ 3x s x b S sb 4tt 

Rosob = Ml = h -T- + — 5 sb 5(x) 

\ r° r 3 3 

Rosab = = RabQs , 

tz z z z ^ f 2M 87r r/^ 
Rpqab = {OpaOqb ~ pb qa ) I — — + — d{X) 

o n/r ( Z X b%q c x b x p c X a Xq X a Xp\ 
-3M [bpa— - \a— - °pb— + Oqb-^) , 

and thus i? o = 4:irG M5 3 (x), R ab = AttG M 5 ab S 3 (x). The solution can be translated and boosted to obtain a moving 
source at an arbitrary location. 

2. Magnetic mass 

We now consider the dual solution, that is, a point magnetic mass sitting at the origin. We have <d 00 (x°,x) = 
N5W(x) as the only non-vanishing component of the magnetic energy-momentum tensor. Furthermore, T^ v = 0. 
The solution is linearized Taub-NUT [20], with only magnetic mass, i.e., in polar coordinates, 

ho<p = —2N(1 — cos#), other components vanish. 

With this choice of ho v the string must be taken along the negative z-axis in order to cancel the singularity at 9 = ir. 
The tensor <& a ^ x is given by $ 0z = — 16nN9(— z)S(x)S(y) (other components vanish). 
One then finds the only non-vanishing components (in Cartesian coordinates) 

Here, Y' a P differs from Y a ^ by a gauge transformation (11.14) with z lm — e lmp ho p , z 0m — 0, and hence gives the same 
curvature. Dealing with Y' a P rather than simplifies the computations. It follows that the curvature is given by 

^OsOb = 0, Rlmab = 0, 

/ 3x b x s S bs 4tt \ 

RlrrMb = Jyeims I ~^ ~ "y bs d(X) 1 , 

U OAT ( 1 , ^X(-\\ Q AT f XbXk XaXk 
ttnmab = ^^abm \ ^3 + — 0\X) I ~ \ €mak ~^5 €mbk 



which satisfies the equations of motion, R a p — and Roijk + RijOk + RjOik = ^TvNeijkS(x) = —8ne ai jx&\- 

Finally, one easily checks that the linearized Riemann tensor of linearized Taub-NUT is indeed dual to the linearized 
Riemann tensor of linearized Schwarschild. In that respect, the reason that it was more convenient to work with Y' a ^ 

instead of Y af ^ above is that it is Y^£ g 7 that is dual to Y^ hw . While the curvatures are dual, the original quantities 
F a/ ^ are dual up to a gauge transformation (11.14). 
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III. MAGNETIC SOURCES FOR BOSONIC HIGHER SPINS 



We now indicate how to couple magnetic sources to spins greater than two. The procedure parallels what we have 
just done for spin 2 but the formulas are somewhat cumbersome because of the extra indices on the fields and the 
extra trace conditions to be taken into account. We shall therefore only give the general outline of the analysis in the 
text and relegate the technical details to appendices. The formalism describing higher spin fields in the absence of 
magnetic sources is recalled in Appendix A. 

The spin s curvature R tllUl n 2 v 2 ---n s v s is the gauge invariant object in terms of which we shall first write the equations 
of the theory. Its index symmetry is described by the Young tableau 



i.e., 



and 



R 



ll\V\~-\XiVi-'ll s V s — R(Jj1V\-- -VilAi" • fast's 1 i — lj * ' * ) S 

^\i\v\--\^iVi^iiJ r {\---^i s v s 0, i = 2, • • • , s — 1 . 



(III.1) 

(III.2) 
(III.3) 



Its dual, defined through 



1 



UP" 



has the same symmetry structure. Note that, just as in the spin 2 case, this does not define an irreducible representation 
of the linear group. But, also as in the spin 2 case, we shall find that only the irreducible part described by 



(III.4) 

(i.e., fulfilling also (III. 3) for i = 1) corresponds to the independent degrees of freedom (the rest being determined by 
the sources). 

The electric and magnetic energy- momentum tensors will be denoted by t fJil ^ 2 ...^ s and /Jil ^ 2 ...^ s . They are conserved, 
i.e. divergence-free: d fl t^ Ul '"' /s - 1 = 0, d tl 9 IJ,Vl "' Va - 1 = 0. Their double traceless parts are written T M1M2 ... Ms and 
MlM2 ... Ms , and are the tensors that actually couple to the spin s field. 



A. Electric and magnetic sources 

The equations in the presence of both electric and magnetic sources are determined again by the requirements: (i) 
that they reduce to the known equations with electric sources only when the magnetic sources are absent, and (ii) 
that they be invariant under the duality transformations that rotate the spin s curvature and its dual, as well as the 
electric and magnetic sources. 

Defining 6^ 1/i2 ... /Js = Q l ^ 1 p, 2 ---p s ~ 1 7 ?(^i^2®^ 3 ... (lts ) > one nn ds the following set of equations for the curvature: 



R n VlV2 

-R[/Jl!^1^2]^2---A t s^s 



d[aRpiv 1 ]ti2i'2---tJ.sV s 



2 -'a»1A»2[A»3[-"[M»i 1/ »]-"] i/ 3] ' 

3 e <wip^ MM-K.^l-H]^] • 



(III.5) 
(III.6) 
(III.7) 



The first equation is the analog of the Einstein equation (II. 5), the second is the analog of the modified cyclic identity 
(II. 1), while the third is the analog of the modified Bianchi identity (II. 3). It follows from these equations that the 
dual curvature obeys similar equations, 



° P\V\H2V2---V.sVs '/ 



">\H\V\IJI,2\v2---llsVs 



9 ®HllJ-2[H3[---[Hs,Vs 



- g € ^lti2P TP V2 [^ 3 [...[^ e , Vs ]...] V3 ]^ 

= g e a^iP^"[^ 3 [.-[n a ,u a ]-]u3]u2]' 



(III.8) 
(III.9) 
(111.10) 



exhibiting manifest duality symmetry. 
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where 



B. Decomposition of the curvature tensor 

As in the spin 2 case, the curvature tensor can be expressed in terms of a completely symmetric potential h^ 1 ... l _ le 
and of a tensor < I )Pcr M1 ... M3 _ 1 fixed by the magnetic energy-momentum tensor, so that the cyclic and Bianchi identities 
do indeed become identities. 

Let $ pcr Ml ... Ms _ 1 be such that 

0,*'V.. /1 .- 1 =*V-m.-i' ( IILn ) 

and let $ p<T AI1 ... AIs _ 1 be the part of $ p<7 ^ 1 ...^ s _ 1 that is traceless in the indices [i\ ■ ■ ■ fi s -i- For computations, it is 
useful to note that d p ^ pcr f _ tl ... f _ tg _ 1 — & (7 ^ 1 ... IXs _ 1 — ^^ r l(mn2®" 7 ^ 3 —^ s _ 1 ) • The expression of the tensor ^..i in 

terms of the Dirac string is given below. The appropriate expression of the curvature tensor in terms of the spin s 
field and the Dirac string contribution is: 

RlllV ll l 2 V2---IJ.sVs = — 2 £ Ml^lpff^ P [n2[H3[—\fla,Ve]—]"3]v2\ ' (III. 12) 

YprT ni-Vs-i = YP<7 H1-H3-1 + ~~ ^' S \pi Y 1 n2-ns-i)e > (III. 13) 

^V--^ = 90X^1...^+^^...^, (111.14) 

* PCT V-^ = - 3(, "^" 2 W ^fcU-o/» • ( IIL15 ) 

The split ofY pcr fll ... fls _ 1 into an A-part and a $-part defines a split of the Riemann tensor analogous to the split (II. 7) 
introduced for spin 2. The Dirac string contribution ($-term) removes the magnetic terms violating the standard 
cyclic and Bianchi identities, leaving one with a tensor r lill/1 ^ 2V2 ... lleVB that fulfills 

and thus derives from a symmetric potential (the spin s field h tll ...^ s ) as 

7 ' P.\V\P.1V1---11 S VB — ' l [jil[f*2"-[f*s,I'a]-"«'2]fl] (III. 16) 

(see Appendix A). The A-term in the curvature is a rewriting of (III. 16) that is convenient for the subsequent analysis. 
The potential h Pl ... fls is determined from the curvature up to a gauge transformation with unconstrained trace. The 
fact that only $ p<T Ml ... Ps _ 1 appears in the curvature and not $ pCT ^ 1 ... / ^ s _ 1 is a hint that only the double traceless part 
••• J u s of the magnetic energy-momentum tensor plays a physical role. 

C. Equations of motion for the spin s field 

In terms of the potential, the remaining equation (III. 5) is of order s. In the sourceless case, one replaces it by 
a second order equation written first by Fronsdal [11]. This can be done also in the presence of both electric and 
magnetic sources by following the procedure described in [14, 15]. The crucial observation is that the curvature is 
related as in (A. 5), namely, 

Rp,\V\^Vl---^sVs r l 1 2 = ~l^y.\tl7\p.3\—\Ha,V a \ — \vz\l (III. 17) 

to the generalized Fronsdal tensor given by 

*Vi-7. = ~\ e ix\xv\ (d X Y»\ 2 ... la -(8-1) a (72 ^ A 73 ... 7s) ) , (111.18) 



so that (III.5) is equivalent to ^^^^[...[^^^...j^j+T^^^f...^^^]...]^] = 0. This implies F^ lfl2l ^ 3 ...^ s +T IJ/llJ/2 ^ 3 ..^ s = 
d(viti2V3^-t*4—t*s) f° r some ^H4---Hs [22]- By making a gauge transformation on the spin s field, one can set the right-hand 
side of this relation equal to zero (see Appendix A), obtaining the field equation 

-F 1 J UlP2M3 - Ms "I" -^f*lf*2A t 3—A»s = 0j (III. 19) 



11 



which fixes the trace of the gauge parameter. When s = 3 this is the end of the story. 

For s > 4 additional restrictions are necessary, namely, we shall demand that the gauge transformation that brings 

the field equation to the form (III. 19) eliminates at the same time the double trace of the field h fll ... IJ/e (see [14] for a 
discussion). 

In terms of the generalized Einstein tensor defined as in (A. 7), i.e., 

G> lM2 ... Ms = F tll p 2 ... lls - ? ?(AiiM2-^) 1 3... Ms )p ' ') (III. 20) 

the equations become 

G> lAI2 p 3 ...p s + 7)ii/J2M.v--Ms = 0. (III. 21) 

We shall thus adopt (III. 21), with the Einstein tensor, Fronsdal tensor and F-tensor defined as in (III. 20), (III. 18) 
and (III. 13), respectively, as the equations of motion for a double traceless spin s field h lll ... lls . These equations 
imply (III. 5) through (III. 10) and define the theory in the presence of both electric and magnetic sources. It is these 
equations that we shall derive from a variational principle. 

D. Point particles sources - Dirac string 

For point sources, the tensors that couple to the spin s field read 

and 

0M"l-"s-l = TVw" 1 



= M/ 1 • • -u^- 1 [ dXS^(x - z(A))i" = M ttVl " n l !— 5^(x-z(x )) 
J u l) 

dXS^ (x - z(X))z" = N - S (3) (x - z{x )) . 



One can check that they are indeed conserved. 

A tensor $ a/3 7l ... 7s _ 1 that satisfies (III. 11) can again be constructed by attaching a Dirac string y M (A,cr) to the 
magnetic source, y M (A, 0) = z M (A). One has 

* Q V-7.-i = Nv ^ ■ ■ ■ / d\da(y"*yP - y a V )S (4) (x - y(X, a)) . 

One can compute explicitly the conserved charges associated with asymptotic symmetries for electric point sources 
(see Appendix A). Using the equations of motion, they read 

where (u) is the traceless part of u Ml • • ■u IJ ">- 1 . One then finds 

/""'""-H") r]^ 1 " 2 ...qfw-if2i u wi+i ... u f.-i)| u |M } 

i 

where the sum goes over I = 0, 1, • • • such that 21 < s — 1, ao = 1 and ai + i = — ^^^gZiI^ a i ■ 
The dual magnetic charges 

are also conserved. 

E. Variational Principle 

The second-order equations of motion G 7l ... 7s + T 7l ... 7s = equivalent to (III. 5), are the Euler-Lagrange derivatives 
with respect to h 11 '"^ 3 of the action 

S[h ltl ... li .(x),y»(\,v)] = J d A x (£ + V-m.^ 1 "'"'). ( IIL22 ) 
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where 

f _ ( S ~ jj y _ T^A I "l Iy «2---a 3 -l I ( S ~ ^) 

2 ^-"-'i 2( S -1)" 

I ( s ~ ^) )iaiyi/pa2-a,-i _ ( S ~ ^) /ia 1 ya 2 pTO3-a s -i 

Indeed, one can check that Sh ^-y s — G 7l ... 7a . As in the spin 2 case, the trajectories of the electric and magnetic 
sources are kept fixed, i.e., the sources are not dynamical. The magnetic coupling in the action was obtained by 
introducing the familiar minimal electric coupling in the "parent action" of [16], which contains two potentials, and 
determining what it becomes in the dual formulation. 

The action (III. 22) can be verified to be invariant under the gauge symmetries (A.l) of the spin s field as well 
as under displacements of the Dirac string (accompanied by an appropriate redefinition of h fll ... fls ). This is done in 
Appendix B. 



F. Quantization condition 

As for spin 2, the unobservability of the Dirac string in the quantum theory leads to a quantization condition. 
The computation proceeds exactly as in the spin 2 case and is reproduced in Appendix C. One finds explicitly the 
quantization condition 

MN 

2^" /-n-^-iW/ 71 '" 7 - 1 («) = «. (ni-23) 

In terms of the conserved charges pTi""7s-i associated with the asymptotic symmetries and their dual charges 
Q7i— 7«-i j the quantization condition can be rewritten as 



1 

2^k 



Q^...^ (w)PT 1 -t- 1 («) e Z . (111.24) 



IV. BEYOND THE LINEAR THEORY FOR SPIN TWO 



We have seen that magnetic sources can be introduced for linearized gravity and linearized higher-spin theories, 
and that an appropriate generalization of the Dirac quantization condition on the sources must hold. However in the 
linear theory the treatment is already incomplete since the sources must be external. In the full non linear theory 
even the introduction of external sources is not possible. For spin 2 the difficulty stems from the fact that the source 
must be covariantly conserved and for spins > 2 the formulation of the non linear theory is still incomplete. 

Nevertheless, we shall address the issue of duality in the spin 2 case, for which the non linear theory in the absence of 
sources is the vacuum Einstein theory of gravitation. This is the "electric" formulation. Electric sources are coupled 
through their standard energy-momentum momentum tensor. We do not know whether magnetic sources should 
appear as independent fundamental degrees of freedom (the complete action with these degrees of freedom included is 
unknown and would presumably be non local, as the results of [23] suggest) or whether they should appear as solitons 
somewhat like in Yang-Mills theory [5]. 

Whatever the answer, there are indications that duality invariance and quantization conditions are valid beyond 
the flat space, linear regime studied above. One indication is given by dimensional reduction of the full Einstein 
theory, which reveals the existence of "hidden symmetries" that include duality [24]. Another indication that non 
linear gravity enjoys duality invariance is given by the existence of the Taub-NUT solution [20], which is an exact 
solution of the vacuum Einstein theory describing a gravitational dyon. The quantization condition on the energy of 
a particle moving in the Taub-NUT geometry is a well known result which has been discussed by many authors [25] 
and which can be viewed as a consequence of the existence of closed timelike lines [26] . 

Because of the importance of the Taub-NUT metric, we shall devote the rest of the paper to a further investigation 
of its properties. First, we shall rederive the quantization condition along new lines, from the quantization of the 
angular momentum stored in the gravitational field. Next, we compute the Poincare charges for the Taub-NUT 
metric, which is asymptotically flat at spatial infinity in the sense of [27]. 
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A. The gravitating magnetic pole 

The Taub-NUT metric is given by : 

ds 2 = -V{r)[dt + 2N(k - cos 9) dif>} 2 + V(r)- l dr 2 + (r 2 + N 2 )(d9 2 + sin 2 9 d<t> 2 ) , (IV.l) 

with 

v{r) = l _ 2(N 2 + Mr) = r 2 -2Mr-N 2 
(> (r 2 +N 2 ) r 2 + N 2 ' 1 ' 

where A and M are the magnetic and electric masses as follows from the asymptotic analysis of the metric and our 
discussion of the linear theory. A pure magnetic mass has M = 0. The number k can be changed according to 



by performing a t coordinate transformation 



k^k' = k-a (IV.3) 



t^t' =t + 2Na<t> . (IV.4) 



The metric (IV.l) is singular on the z-axis. This singularity is most easily seen by calculating |Vt| 2 and can be 
interpreted as a singularity of the metric or of the t coordinate. It is known in the literature as the Dirac-Misncr 
string singularity, due to its analogy with the Maxwell case. Its location depends on the value of k: for k = 1 the 
singularity is at 9 = tt and for k = — 1 it is at = 0. For all other values of k both string singularities exist. The 
choice k = makes the North and South poles play a symmetrical role. In his paper [26], Misner showed that the 
singularity is a coordinate singularity and that the metric describes a non-singular manifold provided that the time 
coordinate t is taken to be periodic with period 8irN. 



B. Spatial rotations and quantization condition 



Even though the metric contains a dtd<p term, it is spherically symmetric. However, the rotation group acts on 
spacetime in an unconventional way [26] and the rotation Killing vectors differ from those of flat space by extra 
terms. One can understand the origin of these extra terms by comparing the Taub-NUT solution with the standard 
electromagnetic magnetic monopole solution and recalling that for stationary metrics, the mixed time-space metric 
components are naturally interpreted as the component of an electromagnetic vector potential ("gravitomagnetism"): 
goi ~ Ai. In this spirit, the above metric component goi would correspond to a monopole potential 

A^ = -2A(fc-cos6>) (IV.5) 

with magnetic charge A, as has been observed by many authors. 

To understand the form of the Killing vectors, we start with the generator of rotations around the z-axis. It is 
recalled in Subsection D 1 that in the electromagnetic case, rotations of the electromagnetic potential of a magnetic pole 
must be accompanied by gauge transformations. When the goi metric components are interpreted as the components 
of a vector potential, the electromagnetic gauge transformations lift to diffeomorphisms along the time direction. We 
therefore expect the gauge parameter which accompanies a rotation to lift to a component along d/dt. 

This expectation turns out to be correct. The metric (IV.l) has four Killing vectors, given in [26] for k = 1, which 
we display as 



e 



d_ 

d „ d „ „„ T cos(A d 



£ x = — sin — — cos 4> cot 9 — + 2 Nk cos <fi cot 9 — 2 A 



89 dej) V " si 1161 / dt 

d ( sin</>\ d 

- - sin^cotfl — + 2Afcsin0cot6» - 2A^^ — , 
o9 dip \ sin 9 J at 

d_ 

dt 



,9 d ( sin</>\ d 

^ y = cos 4> — - sin 4> cot 9 — + ( 2 Afc sin <p cot 9 - 2 A 

6 - — -2Afc-. (IV.6) 



These satisfy the commutation relations 

[&,&] = -e„6c&, [&>,&] =0, (IV.7) 
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where a,b,c — x,y,z and e xyz — 1, which constitute the standard su(2)xu(l) Lie algebra. For N = 0, which is the 
Schwarzschild case, the three Killing vectors £ a generate spacelike two-spheres. In the case of Taub-NUT, although 
the algebra is the same, the action of the group on the manifold is, as pointed out by Misner, different. In fact, the £ a 
which satisfy the su(2) algebra now generate the r =constant three-spheres. These three-spheres have a Lorentzian 
metric. 

The three rotation Killing vectors £ a are invariant under reflections with respect to the origin (t — > —t, — > tt — 8, 
4> — > + 7T, N — > — N and k — ► —k). This should be the case because being the generators of rotations they are 
pseudo- vectors. 

One also observes that £ z has a component along d/dt proportional to the magnetic mass. We find this fact very 
appealing, as it leads to a different derivation of the quantization condition as we now show. 

Consider an electric mass following a geodesic in the Taub-NUT spacetime, with four- momentum p^ = mu^. There 
exists a conserved charge associated with every Killing vector field, and in particular a charge Jz associated with £ z : 

J z = mgun = p$- 2Nmu k . (IV.8) 

We see that just as in the electromagnetic case, the angular momentum about the z-axis has an extra piece besides 
the standard p^. This extra piece comes from the angular momentum in the field, which varies as the particle moves 
and which must be taken into account in the conservation law. In fact, the extra term coincides with the angular 
momentum in the field when k = 1 and the particle is on the positive z-axis (see Appendix D). Thus the angular 
momentum in the field is equal to 2Nmu . Requiring this angular momentum to be quantized in multiples of H/2 
yields the quantization condition 

ANmu Q S KL . (IV.9) 

Note that when the particle is on the negative z-axis, the angular momentum in the field coincides with the extra 
term in Jz when k = — 1 and therefore changes sign. 

Formula (IV.9) agrees with the condition that comes from periodicity in time of the wave function. Suppose indeed 
that the electric mass is described by a wavefunction ip. The time dependence of tp is given by 

ip oc e~ z -yr (IV.10) 

where E = mu is the energy of the particle. Recalling that time is periodic with period At — 8irN, and requiring 
single valucdness of the wave function, we obtain EAt = 2mrh which implies, yet again, that 

ANmu a ehZ. (IV. 11) 

Single-valuedness is required because the closed timelike curves are contractible. 

At this point we may resort to the concept of duality invariance to argue that the result should also hold if we were 
considering a test magnetic mass in the gravitational field of an electric mass, or in fact both. Thus we are led to the 
conclusion that the quantization conditions derived in the linear theory 

h = h G Z (IV ' 12) 

also holds in the full non linear theory. The scalar product which appears here is the Minkowskian scalar product of 
the asymptotic charges Q® and Q bl . The antisymmetry of the quantization condition for the exchange of the dyons 
is confirmed by the change of sign of the angular momentum stored in the field which was observed above. 



C. Taub-NUT is asymptotically flat at spatial infinity 

The introduction of the Minkowskian scalar product is legitimate because the Taub-NUT metric (IV. 1) is asymp- 
totically flat at spatial infinity in the sense of [27]. In polar coordinates, the conditions of [27], including their parity 
conditions, read h rr = 0(r~ 1 ) 1 h r g = 0(1), h r< p = 0(1), hgg = 0(r), hg^ = 0(r), = 0(r) and the leading terms 
(of the order just written) in h rr , h r< f,, hgg, should be even under the inversion 6 — > tt — 6 and <ft — * 4> + n i while 
the leading terms in h r g, hg^ should be odd. Here, hij is the deviation from the flat metric also written in polar 
coordinates. Similarly, the conjugate momenta should fulfill 7r rr = 0(1), ir re = 0(r _1 ), Tr r ^ — 0(r _1 ), tt 00 = 0(r~ 2 ), 
tt ^ = 0(r~ 2 ), = 0(r~ 2 ), with the parity conditions that the leading terms in n rr , tt t ^, tt 00 and ir^ are odd, 
while those in tt t0 and tt ^ are even. 
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It is easy to see that all these conditions are fulfilled in the coordinate system with k = 0, where the Taub-NUT 
metric reads 

ds 2 = -V(r)[dt - 2Ncos6d(f>} 2 + V^^dr 2 + (r 2 + N 2 )(d6 2 + sin 2 6 d(f> 2 ) . (IV.13) 

The perturbation hij is not smooth as a tensor on the sphere because of the string singularities at 8 — 0, ir where 
does not vanish, but its behavior as a function of r is nevertherless compatible with the required asymptotic fall-off. 
This guarantees finiteness of all the surface integrals for the Poincare charges because the parity conditions are in 
particular fulfilled. Coordinate systems with k ^ fail to fulfill the parity conditions and for this reason, arc not 
amenable to the treatment of [27]. 

This notion of asymptotic flatness appears to be all that is needed for the action principle and the identification of 
conserved charges. It does not, and need not, eliminate, of course, the geometrical difficulties pointed out in [26] in 
regard, for example, to the probable impossibility of writing the metric at large distances as a smooth perturbation 
of Minkowski space. 

Note that it is not necessary to fix the magnetic mass N at infinity. Competing histories in the variational principle 
may have different N's. This treats M and N on an equal footing, preserving duality [34]. It is at first sight disturbing 
that the Killing vectors (IV. 6) depend on N because it would seem that the asymptotic symmetry group depends on 
the magnetic mass which as we just noted can be varied. There is however no difficulty because the surface integrals 
associated with the N dependent terms in the killing vectors vanish on account of the parity conditions. Indeed as it 
was observed above since N is a pseudo-scalar the angular dependence of the factor that multiplies it in the rotation 
Killing vectors is also odd under inversions as required by the parity conditions. 

If one computes the Poincare charges for the metric (IV.13) using the formulas of [27], one finds 

P° = M, P*=0, Jij=0, J Ol =0. (IV.14) 

The parameter N does not contribute to the "electric" Poincare charges. Observe in particular that the angular 
momentum is zero for the Taub-NUT solution, indicating that the solution may be thought of as describing a dyon 
with the two charges on top of one another. If the magnetic mass and electric mass are separated, there is a net 
angular momentum and the metric belongs to the Kerr-Newman- Taub-NUT [28] family, which describes metrics with 
angular momentum and electromagnetic charges. These metrics fulfill also the boundary conditions of [27]. Other 
asymptotic analysis of the Taub-NUT solution (at null infinity) may be found in [29]. 



V. CONCLUSIONS 



In this paper we have extended to all integer spins the action for magnetic poles introduced by Dirac. The 
construction has two limitations: (i) it deals only with the linearized field theory and (ii) the sources must be 
external. The presumption is however strong that, at least in the spin 2 case, that is, the theory of gravitation, 
monopoles exist and the corresponding Dirac quantization condition holds. Arguments in this direction are advanced 
based on properties of Taub-NUT space. The analogy with electromagnetism is emphasized throughout. 
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APPENDIX A: HIGHER SPIN GAUGE FIELDS IN THE ABSENCE OF MAGNETIC SOURCES 

We recall in this appendix the free theory of bosonic higher spin fields. 

1. Spin s field and gauge invariance 

In the absence of sources, a massless bosonic spin s field is described by a gauge potential which 
is a totally symmetric tensor /i MlM2 ... Ms subject to the "double-tracelessness condition" [11], h il ^ jl2 ... jle = 
h (mtJ,2-ns), ^iP2P3M4---Ma ? 7 AilA12 »7 A13Ai4 = 0- The S au g c transformation reads 

hfj,ni2---n s ~^ h>mn2---ns + ®{h-l^h2---hs) ' (A-l) 

where the gauge parameter £, fl2 ... fls is traceless, i^^---^^ 2 ^ 3 = 0. The trace condition on the gauge parameter 
appears for spins > 3, while the double tracelessness condition on the field appears for spins > 4. 

From the field h ilxji2 ... ile , one can construct a curvature Rn lvl ii 2V2 —n sVs that contains s derivatives of the field and 
that is gauge invariant under the transformations (A.l) even if the gauge parameter is not traceless, 

RniVlH2V2---HeVs = _ 2 ft'[^i[^2-"[;*»,I'«]-"f2]fl] I (A-2) 

where one antisymmetrizes over and Vk for each k. This is the analog of the Riemann tensor of the spin 2 case. 
The curvature Rm^^w-neVs has the symmetry characterized by the Young tableau 



(A.3) 



and fulfills the Bianchi identity 

^[a-R/J,! I^l]/^2^2 " 'Ms v s 0- (A-4) 

Conversely, given a tensor i? Mll / lAt2 ^ 2 ... MsI / s with the Young tableau symmetry (A.3) and fulfilling the Bianchi iden- 
tity (A. 4), there exists a "potential" h IJ/lfl2 ...^ e such that (A. 2) holds. This potential is determined up to a gauge 
transformation (A.l) where the gauge parameter ^ tl2 ... t _ ls is unconstrained [22]. 

2. Equations of motion 

The trace conditions on the gauge parameter for spins > 3 are necessary in order to construct second-order invariants 
- and thus, in particular, gauge invariant second-order equations of motion. One can show that the Fronsdal tensor 

F — 3 P 8 h - h3, 3 p h > 4- S<yS - — 8. h 9 

which contains only second derivatives of the potential, transforms under a gauge transformation (A.l) into the trace 
of the gauge parameter 

(*-!)(* -2) 

^MlM2---Ms ~~ ' - r MlM2'"Ms ' 2 ty (/ J iA'2M3S A j 4 ... AIs )p ; 

and is thus gauge invariant when the gauge parameter is requested to be traceless. The Fronsdal tensor is related to 
the curvature by the relation 

Rfl 1 U 1 fl 2 V2---HsV s r l 1 2 = _ 2"^1M2[^3[---[Ms,^s]--']^3] • (A-5) 

In the absence of sources the equations of motion that follow from the variational principle written in the main text 
(111.22) are 

<W..„. =0, (A.6) 

where the "Einstein" tensor is defined as 

Gmn2---na — ^ 7 miM2---ms ^ r l(pip, 2 F fl3 ... lls ) p 9 ■ (A-7) 
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They obviously imply 

R^ivni 2 v2---^sv s r l 1 2 = 0, (A. 8) 

and the inverse implication is true as well [15]. Indeed, Eq. (A. 8) implies that the Fronsdal tensor has the form 
Ff ilf _ l2 ... f _ ls = 9( /JlA12M3 E Al4 ... Ats ), which can be made to vanish by a gauge transformation with an unconstrained gauge 
parameter (see [14] for a discussion of the subtleties associated with the double tracelessness of the spin s field 
hfj, 1 ... f j, s ). The interest of the equations (A. 8) derived from the Einstein equations is that they contain the same 
number of derivatives as the curvature. Thus, they are useful to exhibit duality, which rotates the equations of motion 
and the cyclic identities on the curvature. 



3. Dual curvature 



The dual of the curvature tensor is defined by 

and, as a consequence of the equations of motion (A. 8), of the symmetry of the curvature and of the Bianchi 
identity (A. 4), it has the same symmetry as the curvature and fulfills the equations S illVimv ^... ilsVs Tf lV ^ = 0, 

9[ a ^fi 1 u 1 ]fi 2 ' J 2---tJsVs 0- 



4. Conserved charges 

Non vanishing conserved charges can be associated with the gauge transformations (A.l) that tend to Killing tensors 
at infinity ("improper gauge transformations"). They can be computed from the constraints [27] or equivalently from 
the knowledge of their associated conserved antisymmetric tensors k^"^ introduced in [30], which generalize the 
electromagnetic Fp V , the divergence of which vanishes in the absence of sources. The corresponding charge is given by 
= i [ Q /3] dx a A dx 13 , where the integral is taken at constant time, over the 2-sphere at infinity. The tensors 

k^"® read 

u\- a P} — f}a h PtJ 1 ---tJ.s-l p ( S ~ jj a/3}, PHi---fJs-2Ca 

K (, ~On + 2 P ^ A»l-f*a-2 

+(s- l)d p h pa ^-^- 2 ^ u - ^ ~ ^ d^h^-P-'-^P^P u , 
-(a<-/3) + --- , 



where the dots stand for terms involving derivatives of the gauge parameters. 

Of particular interest are the charges corresponding to gauge transformations that are "asymptotic translations" , 
i.e. — ►j.^oo e Mi - Ma-i f or some traceless constant tensor e vi---Hs-i _ p or these transformations, the charges 

become, using Stokes' theorem and the explicit expression for k^"® , 

Q, - I ' " >l'\r. 

Jv 

As these charges are conserved for any traceless e lil ... tie _ 1 , the quantities pA*i"-A*«-i defined as the traceless parts of 
J v G°P' 1 '" fJ - B - 1 d 3 x are conserved as well. In the spin 2 case, P M is the energy-momentum 4- vector. 

5. Electric sources 

In the presence of only electric sources, the equations of motion read: 

Gfiifj 2 ---fJs "b r Ff_ L1 f i2 ...f ls =0, (A. 9) 

or equivalently 

RfJ 1 V 1 fJ2V 2 ---[JsVs r l — 2 r ^fJlfJ2[fJ3[---[fJs,V s ]---]v3] (A. 10) 
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where T, 



T, 



-vi T' 



generalizes the energy momentum tensor and primes denote traces, 
T'^ 3 ... fls — T^...^^ 1 ^ 2 . The curvature tensor has the Young symmetry (A. 3) and fulfills the Bianchi identity (A. 4), 
as in the case without sources. 

On the other hand, while the trace of the dual curvature tensor still vanishes, the latter has no longer the Young 
symmetry (A. 3) and its Bianchi identity gets modified as well. The new symmetry is described by the Young tableau 



I'- 2 






"2 




v s 



as the dual curvature now satisfies S\ 
becomes dr^S 



If T P 



[m»>"«] 



■[fi. 



(A.11) 

while the Bianchi identity 



■]"3]f2] 



APPENDIX B: GAUGE INVARIANCES OF HIGHER SPIN ACTION (WITH MAGNETIC TERMS) 

We prove in this appendix that the action (III. 22) is invariant both under the standard symmetries of the spin s 

gauge field and under displacements of the Dirac string. 

To that end, we first observe that the first term in the action (III. 22) is invariant under the following shifts of 
yP v 



Ql—ae-l ' 



^ Q ,.. as _ 1 = W-.-X) - ^ W-^o + <*«i*"V.«._ l) ■ ( R1 ) 

where z^ v ai ... ae _ 2 — z Lllv \ 1 ... as _ 2 — z^" \ ai ... a ^ is an arbitrary traceless tensor that satisfies rj ai ^ x z IJ '^ ai ... aa _ 2 = 

when s > 2. Under this transformation, Y^ v Ql ... Q . 3 _ 1 transforms as SY^ V ai ... ae l — d^z 111 ' a2 ... a t y which makes it 

obvious that the curvature and the Fronsdal tensor are invariant under (B.l). 
The transformation (B.l) can be conveniently written 

^ ai ... at _, = e^dpa^...^ , (B.2) 
where a aai a 2 - as _ 1 = -a aiaa2 ... ole _ 1 = o CTai(a2 ... Qs _ 1 ) is given by 

1uaia2-a 8 -i = ^ e f/37"i z ^ ' ol2 ... ote _ 1 , (B-3) 

is traceless and satisfies a[ CTQlQ2 ] Q3 ... Q3 _ 1 = when s > 2. 

Direct computation shows that the gauge transformation (A.l) of the spin s field acts on ai ... as _ 1 as the 
transformation (B.2) with 

= ^-^ B P 

(Jj pcr{oLi---cxs-2) g (J [p<,a]ai---a s -2 

i (s - !)(« - 2) 



It follows from this fact and the conservation of the energy-momentum tensor that the action (III. 22) is invariant 
under the standard gauge transformation (A.l) of the spin s field. 

The displacements of the Dirac string change ^ v ai ---a s - 1 as S^ v ai ... cte _ 1 = k^ v ctl ... cte _ 1 where d l ^k^ v ctl ... as _ 1 = 
0. The latter equation implies that k^ v ai ... as _ 1 — d\K tlvX <Xl ... as _ 1 , where K' luX ai ... as _ 1 — K^ vX } Xl ... as _ 1 . Let 
^"oi-a,-! be the part of K^ vX ctl ... as _ 1 that is traceless in a\ ■ ■ - a s -i; it can be decomposed as 

ai-«8-i ai-o e -i ~ (ai y a 2 ---a 3 -i) ' 

where _! satisfy 

^P-vX _ [ftvX] _ uuX uv\ mi _ r, 

x ai-a s _i _ x ai-a,_i— ( ai ... Qs _i)> ai-a s _i°A ~ u > 

f a 2 ---a 5 _i — f 02-q,-i — J/ (aj.-a,.!) i f a 2 ---a.,-i t V ~~ u ' 

— _ (f__jO jlii/A „aia 2 \xv rl a 1 a 2 _ n 

'/ </ ai-a s _i - ai-a,_i'/ )» ni-a,_i'( — U. 
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For the action to be invariant under displacements of the string, the variation of 3> All ' ai ... Qs _ 1 has to be supplemented 
with an appropriate transformation of h ai ... as . This transformation reads 5h ai ... as = ^ e^AOi X ^ V a2 ---a ) ■ Indeed, 
when one performs both variations, Y^ v ai ... cts _ 1 transforms as in (B.l), so the first term in the action is invariant. 
Furthermore, the electric coupling term is invariant as well because the support of the variation of the spin s field 
does not contain the electric worldlines. 

The identities that follow from the invariance (B.l) - or (B.2) - of the first term C in the action may be written 
conveniently in terms of 

^0-71 -73-1 _ e °V-v\Q ( ^ \ m 4 x 

and its trace A'T 2 "'^- 1 = A""* 1 "^'- 1 ^^ . They read 

= A a " <1 '"' 1a - 1 — A^ 1 ' l2 '" la ~ 1 ^ a — S — - (rf^ 11 ^'^'"T 3 - 1 ) _ jy(7l72^'73---7s-l)o-^ 5-j 

Using these identities, one checks the following useful relations: 

Q71--T8 — ^ — ^(7172-"7s) _|_ ( S ~ ~ ^) ^(7172^/73---7») (B.6) 

8h 11 ... ls 2s 

— ^47l72---7s _|_ ( S ~ I)( S ~ ^) ^(7273^74---7s)7l 7) 

which will be used in the following appendix. 



APPENDIX C: QUANTIZATION CONDITION FOR ELECTRIC AND MAGNETIC HIGHER SPIN 

SOURCES 

We work out explicitly in this appendix the quantization condition for higher spins. As for lower spins, this is 
achieved by integrating the constraints expressing the unobservability of the Dirac string. 
In the gauge y° = A, these constraints read 

BC 

7r m = -2Ny' n f 7l ...^_ 1 (v 



u 7i-"7»-i 



In the quantum theory, the wave functional ip must thus fulfill 



7l-"7s-i 



Integrating this equation along a path that encloses an electric source, one finds the following variation of the phase 
of the wave functional 

AT I' BC 

A* = - ¥ / 7l ... 75 _» J ^ ^ y m y' n -y n y' m )dadX, 



where the integral is taken on the two-dimensional surface enclosing the electric source. Using the Gauss theorem, 
this can be converted into a volume integral, 



A* = ~ / 7l ... 7s _» J d 3 x e mn *d p ^ 



BC 



u± 7i-"7s-i 



Using the relations (B.7), one checks that 
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where the dots stand for terms of the form rf 1 ' ll2 X~ <3 '"' 1s - 1 \ Upon use of the Einstein equations G° 71 ■"7»-i = 
_ T o 7 i-7a-i ) the variation of the phase becomes, A* = f / 71 ... 7s _ 1 (v) J d 3 x T^-t- 1 = f- fl ...- /e _ 1 ( v )F 1 "" r '- 1 ( u )- 
For the wave functional to be single-valued, this should be a multiple of 27T. This yields the quantization condition 

MN 

2^-/7i-7.-i(«)/ 71 '" 7 - 1 («) = «. (C- 1 ) 
Introducing the conserved charges i J 7 1- "7»-i ) Q71— 7«-i j this can be rewritten as 

lo 1r .. Vl (^'''^( u )eZ. (C.2) 

APPENDIX D: ROTATION SYMMETRY AND CHARGE QUANTIZATION FOR AN 
ELECTROMAGNETIC MAGNETIC POLE 

1. The Dirac string and rotations 

A monopolc is spherically symmetric, yet the vector potential is not manifestly so, and when rotated, the string 
will move to a new location. As is well known, spherical symmetry is recovered by supplementing the naive rotation 
generated by the usual Killing vectors by a gauge transformation that moves the string back to its original location. 
Let us see how this works explicitly. 

The usual Killing vectors which generate SO(3) rotations on the two-sphere are given by 

6f = — sin — - cos0cot6» — , 

. d . 1 d 

£y = cos — - sin <p cot — , 
a0 dip 



They satisfy the algebra 



[U,&] = -e A Bctic- (D.2) 



Under an infinitesimal rotation generated by one of these vectors, the vector potential changes by the Lie derivative 
along the vector field: 

C,,A, = ZidjAi + AMb = e B Foi + • (D.3) 

For the monopole vector potential (D.ll) this Lie derivative is non-zero and hence usual rotational symmetry would 
appear to be broken. However, we can undo its effect if we accompany the Lie derivative by a gauge transformation 
Ai — > Ai + diA such that the total change is 

5 iB A i = C SB A i +d i A B = 0. (D.4) 
The gauge transformations that we find associated with the three Killing vector fields are 

9 _„^„ : _ fl c j a 9k_, n _„ ± 5 cosi 



Ax = — — cos</>sin# — £, x Aj = — cot 6 cos 



47r A J 47T 47r sin6» 

9 j.^ a e j a - 9 k jl 9 sin< 



Ay = — — sin sin — ZyAj — — cot 8 sin< 



47T Y J 47T 47rsin0 

^— cos 8 — £tAi = — ^— 
4tt sz 3 4tt 



A 2 = -TZ<x»0-& A i = -7Z- ( D - 5 ) 



We note that Az is a constant gauge parameter, yet it needs to be considered in order that the new rotation generators 
to be defined below satisfy the algebra. 

We can combine the gauge transformations A# and the naive rotation operators £b to create new rotation operators 
£b = Cb + Ab^Jj by using the fact that the generator of gauge transformations is the electric charge q, which is in 
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turn conjugate to the coordinate on the U(l) fiber which we call A, so we can write: 

" . , d n d ( gk . g cos<j>\ d 

E x = — sin 6 — — cos cot — — + — cos <b cot — 



de Y d0 1 V4tt y A^smejdx 1 

. gk . , g sin</>\ 9 

4y = cos — — sin <p cot & — + — sin cot ( 



COS' 




« 9 
6cot6> — 

00 


d 


gk d 




~d4> 


4tt d\ ' 





At: y 47rsin0yaA' 

** = -B$-%tx- (D ' 6) 

These operators satisfy the SO(3) algebra just as (D.l). 

2. The field angular momentum 

Let us now consider an electric charge q of mass m placed at a point a in the background field of a magnetic pole 
g placed at the origin. The background magnetic and dynamical electric fields are 

Q x q(x — a ) . 

47r|x| 3 Att\x - a\ 3 y ' 

As is well known, the angular momentum stored in the electromagnetic field does not vanish. It can be calculated 
from the standard symmetric energy-momentum tensor of the electromagnetic field and is given by 

Lf ie i d = [ xx(ExB) = ^L£- . (D.8) 
Jv 47r l a l 

This expression is independent of the magnitude of a, though it depends on its direction, and therefore the limit 
a —> cannot be taken continuously. For example, if we place the charge above the pole the angular momentum will 
be negative, whereas it will be positive if placed below. Whatever the sign, if we assume the angular momentum to 
be quantized in integer multiples of h/2, we obtain the quantization condition 

^£Z. (D.9) 

This was first noted by Fierz and is the statement that quantization of the field angular momentum results in the 
Dirac quantization condition. The fact that one may obtain half integer angular momentum without fermions in sight 
has been lucidly analyzed in [31]. 

3. Angular momentum of a test electric charge in the background of a magnetic monopole 

One can get direct access to the angular momentum in the field differently by analyzing the angular momentum of 
a test electric charge in the background of a magnetic pole. 
The magnetic field of a magnetic monopole of strength g, 

r-JL,, (D.io) 

can be associated with the one form vector potential 

A= JL (k- cos 8) d<p, (D.ll) 

which satisfies B = V x A, for k an arbitrary constant, away from the Dirac string. When k = 1 the Dirac string lies 
along the negative z-axis while when k = — 1 the Dirac string lies along the positive z-axis. The value of k can be 
modified to k by a gauge transformation with parameter g(k — k)4>/4ir. For values of k ^= ±1, there are two strings, 
one along the positive z-axis and one along the negative z-axis, with strengths that add up to one. 

By demanding that the velocity should transform as a vector under rotations, the angular momentum operator 



(D.12) 
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was derived in [32] . 

There is an alternative way to arrive at this result which we believe is worth recalling. Using the standard Noether 
procedure for the action 



/G 



^miiX 1 - qAiX 1 j dt (D.13) 

we obtain that for a variation Sx the action changes by 

5S = j j t [(mi, - qA t )Sx l ] . (BAA) 

on-shell, i.e. when the equations of motion hold. The transformation will generate a conserved charge if the action is 
invariant off-shell up to the integral of a total derivative. Let us consider what the case is for a rotation generated by 
Sx 1 = £g. The first term in (D.13) is the norm of a vector and therefore manifestly invariant under rotations. The 
variation of the second term becomes 

SS = -q j (x i 6 iB A i + A.S^x 1 ) 

= -qj {C iB Ai)x\ (D.15) 

where we have used S^ B A{ — (djAi) 5^ B x^ . If, as is the case here (see (D.4)), the Lie derivative of A t is given by a 
total derivative C^ B A t = -diA B this becomes 



5S = qjj t [K B \. (D.16) 

We can now equate equations (D.14) and (D.16) to obtain the charge J B conserved under rotations 

J B = (mxi-qA i )t B -qA B - (D.17) 

In fact the A B are the gauge transformations we computed in (D.5). For example, for rotations about the z-axis, we 
obtain that 

J z = m<j)+ ^cos(9, (D.18) 

and in general the result (D.12). 

From the expression (D.12), it is clear that the angular momentum docs not vanish when the particle is at rest. 
This apparent paradox is explained by comparing (D.12) with (D.8). The extra piece is the angular momentum of the 
electromagnetic field, which must be included because the electric field of the test particle is dynamical and hence, 
the angular momentum of the electromagnetic field changes with time. Only the sum of the standard orbital angular 
momentum and the field angular momentum is conserved. 

Expressing the angular momentum in terms of the (non gauge invariant) conjugate momenta p<f, = m<j) — qA^ one 
obtains, for the angular momentum about the z-axis, 

Jz=Pcb + k^-. (D.19) 

In the gauge k = 1 and when the particle is on the positive z-axis, the field angular momentum is just equal to the 
extra piece appearing in Jz, besides the usual p^. This is a way to identify the field angular momentum if one knows 
the total angular momentum, as it is the case for gravity. Note that in the gauge k — 0, p^ is equal to the total 
angular momentum. 

It has been shown in [31, 32], that the difference Jz — qg/^ has integer eigenvalues. We may rewrite this assertion 
in the notation used here as the statement that in the gauge k = 1 the operator P4> — j-§^ has integer eigenvalues 
and therefore the wave function is periodic in <j). 

However the periodicity of the wave function depends on k. This can be traced to the fact that differentiating with 
respect to (f> 
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at constant A (where A is the coordinate along the U(l) fiber) is not invariant under the gauge transformation 




(D.21) 



which brings k from 1 to 0. As (D.21) itself shows the wave function picks up a phase factor e~ 1 ^^ which shows in 
turn, that in the gauge k = the wave function is antiperiodic in <j> when qg/Air is half integer. This anti-periodicity 
is permissible because in the gauge k = an infinite line is removed from R 3 and the resulting configuration space for 
the electron is therefore not simply connected. The infinite line is formed by two strings of half strength that come 
into the magnetic pole along the positive and negative sides of the z-axis. 

There is another interesting manifestation of the lack of simple connectedness due to the infinite line formed by the 
two strings when k = 0. It is the following: the closed path traveled in configuration space when the string sweeps 
out a closed surface around the electric charge is not contractible to the identity, whereas it is so when the string is 
half infinite (see [10]). For this reason the string wave function changes sign after the turn and one gets the same 
quantization condition as when k = 1 in which case the path is contractible and the wave function returns to its 
original value. 
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